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The time derivative of the kilohertz quasi-periodic 
oscillations in 4U 1636—53 
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ABSTRACT 

We analysed all archival RXTE observations of the neutron-star low-mass X-ray binary 
4U 1636-53 up to May 2010. In 528 out of 1280 observations we detected kilohertz 
quasi-periodic oscillations (kHz QPOs), with ~ 65% of these detections corresponding 
to the so-called lower kHz QPO. Using this QPO we measured, for the first time, the 
rate at which the QPO frequency changes as a function of QPO frequency. For this we 
used the spread of the QPO frequency over groups of 10 consecutive measurements, 
sampling timescales between 320 and 1600 s, and the time derivative of the QPO 
frequency, i>qpo, over timescales of 32 to 160 s. We found that: (i) Both the QPO- 
frequency spread and t'Qpo decrease by a factor ~ 3 as the QPO frequency increases, 
(ii) The average value of i'qpo decreases by a factor of ~ 2 as the timescale over which 
the derivative is measured increases from less than 64 s to 160 s. (iii) The relation 
between the absolute value of ;>qpo and the QPO frequency is consistent with being 
the same both for the positive and negative QPO-frequency derivative. We show that, 
if either the lower or the upper kHz QPO reflects the Keplerian frequency at the inner 
edge of the accretion disc, these results support a scenario in which the inner part 
of the accretion disc is truncated at a radius that is set by the combined effect of 
viscosity and radiation drag. 
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1 INTRODUCTION 

Kilohertz quasi-periodic oscillations (kHz QPOs) are the 
fastest variability so far observed in neutron star (NS) low- 
mass X-ray bin ary (LMXB) systems. The s e oscillations were 
first detected l|van der Klis et all Il996al ; IStrohmaver et all 
1 1996T ) shortly after the launch of the Rossi X-ray Timing 
Explorer (hereafter RXT E; see lBradt. Rothschild fc Swankl 
ll993l : |Jahoda et alj|2006h . Since then, kHz QPOs have been 
observed in more than 30 neutron star LMXBs (see, e.g., 
Ivan der Ku1l2005l 1. KHz QPOs are often observed in pairs, 
usually called lower and upper kHz QPO, with frequencies 
rangi ng from a few hun dred Hz to more than 1 kilohertz 
fe.g.. IvarTd cr Klis 20051. for an extensive review). The fact 
that the timescale of these oscillations corresponds to the 
dynamical timescale very close to the NS, makes kHz QPOs 
potentially one of the few tools nowadays available to di- 
rectly measure strong-field gravitational effects. 

In the past 16 years, several models have been 
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proposed to explain these oscillations, with empha- 
sis in trying to explain the frequency of the QPOs. 
Almost all models associate the frequency of these 
oscillations with characteristic fre quencies in a geomet- 
ricall y thin accretion dis c (e.g., iMiller. Lamb fc Psaltid 
1 19981 ; IStella fc Vetril Il998f ). Systematic studies of other 
properties of kHz QPOs, such as the coherence and 

the amplitude of the oscillati on, are also available for 

sever a l LMXBs (IJonker et all l2000l: Ivan Straaten et al 



|2000l ; |Pi Salvo et alj l200ll ; iMendez. van der Klis fc Ford 
' 20011 ; iHoman et all |2002|; Ivan Straaten et al 



2002 



Di Salvo, Mendez and van der Klisl 



2003 



van Straaten et al.l 2003; B arret et al. l2005al; I 3arret et al.l 
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2005bl;lBarret. Olive fc Milleijl2006l ; lAltamirano et al.ll2008l ; 

2010h . but not many models explain thes e 



other properties fsee lMiller et al.lll99Sl ; iBarret et al.ll2006t ). 
Some sources show a characteristic dependence of the 
coherence (defined as Q = v/Av, where v and Av are the 
centroid frequency and the full width at half-maximum of 
the kHz QPO, respectively) and the amplitude of those 
oscillations upon the centroid frequency of the oscillation: 
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Both the coherence and the amplitude of the lower kHz 
QPO in these sources increase as the frequency increases, 
and then q uickly decreas e as th e QPO frequency continues 
to increase. I Barret et alj l|2006l ) proposed that the drop of 
the coherence Q at high frequencies is a direct consequence 
of the disc approaching the innermost circular stable orbit 
(ISCO), where the inner edge of the disc gets shattered 
by the gravitational force of the neutron star. If this 
interpretation is correct, the drop of the coherence and 
rms amplitude of the lower kHz QPO at high frequencies 
would confirm the existence of the ISCO, and it would 
provi de a direct measu rement of the mass of the neutron 
star (I Miller et alj [l998) . However, other non space-time- 
related interpretations h ave been prop osed to explain the 
properties of QPOs; e.g., iMendej (2006) suggested that the 
drop of Q and rms amplitude in individual sources might 
be related (at least in part) to changes in the properties 
of the accretion flow in these systems. This could explain 
the significantly different properties of the kHz QPOs, 
respectively, in the high and low-lumi nosity phase of th e 
outburst of the LMXB XTE J1701-462 ijSanna et al.ll2010h . 

Here we investigate the properties of the kHz QPOs for 
the NS system 4U 1636-53 by scanning the whole RXTE 
archive until May 2010. In Section 2 we describe the obser- 
vations and the data analysis, and in Section 3 we present 
our results. In Section 4 we discuss our results in the context 
of current ideas concerning the origin of the kHz QPOs in 
LMXBs. 



2 OBSERVATIONS AND DATA ANALYSIS 

We analysed all archival observations of the LMXB 4U 1636- 
53 obtained with the Proportional Counter Array (PCA) on 
board of RXTE up to May 2010. We used a total of 1280 
RXTE observations, which correspond to an exposure time 
of ~ 3.5 Ms. During these observations the source showed a 
large number of t ype-I X-ray bursts that we excluded from 
our analysis (see. lzhang. Mendez fc Altamiranoll201lf) . 

Except for four observations, we used event-mode data 
with 125/xs time resolution and 64 channels covering the 
full PCA energy band to produce Fourier power density 
spectra (PDS), setting the time resolution to 4096 points 
per second, corresponding to a Nyquist frequency of 2048 
Hz. In the four observations in which the event-mode data 
were not available, we used a combination of binned modes 
that allowed us to reach the same Nyquist frequency in 
the PDS. In all cases we selected photons with energies 
below ~46 keV (absolute channel 101), with the exception 
of observation 10072-03-01-00 where we selected photons 
between channels 24 and 139, and observations 10072-03- 
02-00 and 10088-01-08-010 where we used the full energy 
band, because in these cases it was not possible to have 
simultaneously the same Nyquist frequency and energy 
selection as in the other observations. We created Leahy- 
normalised PDS for each 16-seconds data segment. We 
removed detector drop-outs, but no background subtraction 
or dead-time correction were applied before calculating the 
PDS. 

We averaged all 16-s spectra within each observation, and 
we searched for kHz QPOs by fitting the averaged PDS in 
the frequency range 200-1500 Hz using a constant to model 



the Poisson noise and one or two Lorentzians to model the 
QPOs. From each fit, we estimated the significance of the 
Lorentzian by dividing its normalisation by the negative 
la error. We considered only features with this ratio larger 
then 3, and coherence Q larger than 2. 

We calculated the X-ray colours and intensity of the source 
using the Standard 2 data (16-s time resolution and 129 
channels covering the entire PCA energy band). We defined 
soft and hard colour as the count rate ratio in the energy 
bands 3.5-6.0 keV / 2.0-3.5 keV and 9.7-16.0 keV / 6.0-9.7 
keV, respectively. The intensity of the source is defined 
as the count rate in the energy band 2.0-16.0 keV. To 
obtain the exact count rate in each of these energy bands 
we linearly interpolated in channel space, since the energy 
boundaries of each channel change slightly with time. To 
correct for gain changes and differences in the effective area 
between the proportional counter units (PCUs) as well as 
differences due to changes in the channel to energy conver- 
sion of the PCUs as a function of time, we normalised the 
colours and intensity to those of the C rab Nebula obtained 
clos e in time to our observ ations (see iKuulkers et al. I ll994l 
and lAltamirano et al.ll2008l for details). Finally we averaged 
the normalised colours and intensities per PCU for the full 
observation using all available PCUs. 



2.1 KHz QPOs Identification 

We detected kHz QPOs in 528 out of 1280 observations anal- 
ysed. Among the observed QPOs, ~84% where detected 
with a significance larger then 3.5a, ~72% with a signifi- 
cance larger then 4ct, and ~53% with a significance larger 
then 5a. Even though ~16% of the QPOs could be con- 
sidered only marginally significant (between 3 and 3.5a), 
all these QPOs follow the frequency- hardness correlation 
described below, which adds confidence to the detections. 
Only 26 observations with QPOs showed two simultaneous 
high-frequency oscillations. In those cases the QPO identi- 
fication is trivial because lower and upper kHz QPOs were 
both present. For all the other observations where we only 
detected a single QPO peak we had to apply a different 
method to identify the peak as either the lower or the up- 
per kHz QPO. In Figure[T]we plot, for each observation, the 
centroid frequency of all the kHz QPOs as a function of the 
hard colour of the source. From the plot it is apparent that, 
for frequencies below 1000 Hz, the data follow two separate 
branches. The first branch extends between QPO frequencies 
of ~ 550 Hz and ~ 950 Hz, with more or less constant hard 
colour at around 0.65, except for a small deviation at low fre- 
quencies where the hard colour increases slightly. The second 
branch shows a clear anti-correlation between frequency and 
hard colour. This has b een already reported for this source 
(see iBelloni et all 120071') as well as for other neutron-star 
syste ms (e.g., IMendez et al.ll 19991 : IMendez and van der Khsl 
1999) and can be used to identify the two different kHz 
QPOs: Lower kHz QPO at low hard colour, and upper kHz 
QPO at high hard colour. 

In FigureQ]we also show the kHz QPOs detected simul- 
taneously per observation: Red-empty triangles and blue- 
empty circles represent, respectively, the lower and the up- 
per kHz QPOs. Those points confirm the correlation previ- 
ously mentioned between frequency, hard colour and type of 
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Figure 1. Centroid frequency of the kHz QPOs detected in 4U 
1636-53 as a function of hard colour. Grey bullets represent the 
QPOs in observations in which only one kHz QPO is observed. 
Red triangles and blue circles represent QPOs identified, respec- 
tively, as the lower and the upper kHz QPO in observations in 
which two simultaneous kHz QPOs are detected. For clarity, in 
the inset we plot only the QPO frequencies of the observations in 
which two simultaneous kHz QPOs were detected. 



QPO. Interestingly, by checking the distribution of the blue 
circles in Figure[T]we can resolve the ambiguity regarding the 
identification of the QPOs in the region above 1000 Hz, in 
which the two branches approach each other. To summarise, 
we detected a total of 357 lower kHz QPOs a nd 197 upper 
kHz QPOs. Compared to iBelloni et atl 120071 we find that 
the branch corresponding to the upper kHz QPO extends 
to higher frequencies at more or less constant hard colour, 
while the branch corresponding to the lower kHz QPO ex- 
tends to lower frequencies at higher hard colour. 



2.2 Variations of the kHz QPO frequency 

It is well known that the frequency of the kHz QPOs can 
change over t ens of Hz in time i ntervals of a few hundred 
seconds (e.g.. iBerger et aTlll996t ). and this can artificially 
broaden the QPO in the averaged power spectrum of long 
observations. In order to study frequency variations of the 
kHz QPOs, we tried to detect kHz QPOs on timescales 
shorter than the full observation. We found that in obser- 
vations with only the upper kHz QPO we cannot detect 
the QPO significantly on timescales shorter than the full 
observation, whereas we detected the lower kHz QPO on 
timescales between 32 s and 160 s in about 50% of the ob- 
servations that showed kHz QPOs, all at significance levels 
larger than 4u. For these observations w e created dynamica l 
power density spectra (see, e.g., Fig. 2 in lBerger et al.| [l996). 
and we extracted the frequency evolution as a function of 
time (hereafter frequency profile). To reduce the influence of 
the frequency error on the analysis of th e frequency varia- 
tion, we applied a Savitzky-Golay filter (|Savitzkv fc Golavl 



1 1964 1 Press fc Vetterlindl 19891 ) to smooth the frequency pro- 
file. This method locally fits a polynomial function on a se- 
ries of values (in our case we used a 4th order polynomial on 
63 consecutive points), and it provides the first time deriva- 
tive of the QPO frequency (hereafter frequency derivative) 
for each measurement. 

For each observation we then calculated how much the 
centroid frequency moved with time by measuring the sam- 
ple standard deviation (hereafter the spread) using groups 
of 10 consecutive measurements of the frequency profile. We 
also used the frequency derivative to evaluate the rate at 
which the centroid frequency changes on certain timescales. 



3 RESULTS 

In Figure [2] we show the spread of the lower kHz QPO fre- 
quency as a function of the QPO frequency measured over 
timescales of 320 and 1600 s. Each point represents the av- 
erage spread within frequency bins of 5 to 40 Hz, depending 
on the number of measurements at different frequencies. No- 
tice that the rebin has been applied under the assumption 
that the spread only depends on frequency. 

Figure[2]shows that the frequency spread decreases from 
~ 7 Hz to ~ 2 Hz as the QPO frequency increases from ~ 
650 Hz up to ~ 910 Hz. Above ~ 910 Hz, the data show a 
marginal increase of the spread to ~ 4 Hz as the frequency 
increases further up to ~ 930 Hz. The apparent rise of the 
spread at high frequency is suggested by two points only, 
which are the average of 13 and 12 measurements, respec- 
tively. We fitted the data both with a line and a broken 
line function, and we carried out an F-test to compare both 
fits. The F-test probability is 2.2 x 10 -2 , indicating that the 
fit with a broken line is marginally better than a fit with 
line. We cannot rule out the possibility that this increase is 
real, but as we discuss below, it could also be (partly) due 
to a larger contribution of the statistical errors to the QPO 
frequency in that frequency range. 

To take into account the contribution of the statisti- 
cal errors from the measurements of the frequency on the 
spread, we subtracted in quadrature the 1-sigm a errors from 
the observed spread (see IVaughan et al.ll2003l , section 6.1, 
Eq. 8) . We found that the errors of the QPO frequency con- 
tribute about 10% of the observed spread at around 650 
Hz, whereas this contribution decreases down to 5% at 850 
Hz, and then increases up to 40% at ~ 930 Hz. Although 
this correction slightly flattens the relation, the overall trend 
does not change significantly, except that the apparent in- 
crease of the spread at the high-frequency end becomes less 
significant. 

The accuracy with which one measures the centroid fre- 
quency of the QPO depe nds on the amplitu de and the width 
of the oscillation (e.g., Ivan der Klidll997h. which in turn 
change with QPO frequency (e.g., Barret et alj|2006h . In or- 
der to verify whether this affects the relation between spread 
and frequency, we made simulations where we added nor- 
mally distributed noise to the frequency profile, with stan- 
dard deviation ranging from 1 to 10 Hz. As expected, the 



1 We varied the number of overlapping points from 4 to 10, but 
within errors this does not affect the time derivative of the QPO 
frequency. 
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Figure 2. Measurements of the frequency spread of the lower kHz 
QPO in 4U 1636-53 as a function of the QPO centroid frequency. 
Each point is the average of the QPO-frequency spread at that 
frequency, corresponding to timescales between 320 s and 1600 s. 
See section [3] for more details. 



Figure 3. Absolute value of the frequency derivative of the lower 
kHz QPO frequency in 4U 1636-53 as a function of the QPO 
frequency. Each point is the average value (within a frequency 
bin) of the frequency derivative of the QPO frequency detected 
on timescales of 64 s or less. See section [2.21 for more details. 



extra noise component tends to flatten the relation, but er- 
rors of up to 10 Hz (rms) in the determination of the QPO 
frequency do not significantly change the trend shown in 
Figured 

Figure [3] shows the absolute value of the frequency 
derivative of the lower kHz QPO as a function of the QPO 
frequency. This plot is based on measurements of the QPO 
frequency on timescales of 64 s or less. Each point repre- 
sents the average value of the frequency derivative within 
a frequency bin from 5 to 30 Hz. Similar to what happens 
with the spread, the frequency derivative decreases as the 
frequency of the QPO increases, but unlike Figure [2] in this 
case there is no indication of an increase of the frequency 
derivative above ~ 910 Hz. 

In Figure [4] we compare the absolute value of the fre- 
quency derivative of the lower kHz QPOs detected over three 
different timescales. Black points represent the frequency 
derivative of the QPOs measured on timescales of 64 s or 
less (the same as Figure |3J, red points are for timescales be- 
tween 64 and 160 s, and green points are for timescales of 
160 s. Interestingly, for timescales longer than 64 s the fre- 
quency derivative initially stays more or less constant as the 
QPO frequency increases, and at around 800 Hz it decreases 
following more or less the same trend seen on the shortest 
timescale. We investigated whether the difference between 
trends in Figure U is only caused by the different timescales. 
To do that we took the observations where we detected kHz 
QPOs in segments of data shorter than 64 s (black points in 
Figure we created PDS increasing the length of the data 
segments up to 160 s and we searched again for kHz QPOs. 
We then recalculated the frequency derivative and we com- 
pared it with the green points in Figure[3] We found that the 
two distributions of points are consistent with each other, 
which shows that the difference between the three groups of 



points in Figure [4] is likely due to the different timescales 
over which we measured the QPO frequency. 

Finally, in Figure [5] we show, respectively, the abso- 
lute value of the negative (frequency decrease) and positive 
(frequency increase) frequency derivative for kHz QPOs de- 
tected on timescales of up to 64 s. We fitted both relations 
with a linear function. We found that, within the errors, 
there is no significant difference between the trend of the 
positive and negative QPO-frequency derivative with QPO 
frequency. 



4 DISCUSSION 

We detected kHz QPOs in 528 out of 1280 RXTE observa- 
tions of the LMXB 4U 1636-53 (all RXTE data available of 
this source up to May 2010). 26 of these 528 observations 
showed simultaneously the lower and the upper kHz QPO, 
whereas in the remaining 502 observations we detected a 
single kHz QPO that we identified using the frequency vs. 
hard colour diagram (see section [2~T)l . Using measurements 
of the lower kHz QPO on timescales between 32 and 160 s 
we measured, for the first time, the rate at which the QPO 
frequency changes as a function of the QPO frequency itself. 
For this we used both the spread of the QPO frequency over 
groups of 10 consecutive measurements, sampling timescales 
between 320 and 1600 s, and the time derivative of the QPO 
frequency over timescales of 32 to 160 s. Both the QPO- 
frequency spread and QPO-frequency derivative decrease by 
a factor of ~ 3 as the QPO centroid frequency increases from 
~ 600 to ~ 900 Hz. Above ~ 910 Hz the QPO-frequency 
spread appears to increase again, albeit the evidence for 
this is marginal. We found that the relation between QPO 
frequency, i^qpo, and QPO-frequency derivative, vqpo, de- 
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Figure 4. Absolute frequency derivative values of the lower kHz 
QPO in 4U 1636-53 for 3 different timescales. Black points are 
measurements for QPOs detected on intervals shorter than 64 s. 
Red points represent detections between 64 s and 160 s, while 
green points represent measurements of the QPO on intervals of 
160 seconds. 



pends upon the timescale over which the QPO frequency is 
measured: When ^qpo is between 600 and ~ 770 Hz, z>qpo 
decreases by a factor of ~ 2 as the timescale over which we 
measure the frequency derivative increases from 64 s or less 
up to 160 s. Finally, we found that the positive and negative 
values of vqpo are consistently distributed in the whole fre- 
quency range (see Figure [5]), even at high frequencies where 
the accretion disc shoul d approach the ISCO. 

IBarret et ail (|2006t ) found that in 4U 1636-53, the co- 
herence of the lower kHz QPO increases from ~10 to ~220 as 
the frequency of the QPO increases from ~600 to ~850 Hz, 
and then it decreases rather abru ptly as the frequenc y con- 
tinues increasing up to 950 Hfl IBarret et alJ lj2006T ) were 
able to model the behaviour of the coherence of the kHz 
QPO in terms of the lifetime of blobs of matter that move in 
quasi-Keplerian orbits at the inner edge of a geometrically- 
thin disc, the width of the disc annulus where these blobs 
produce the kHz QPOs, and the advection speed of the gas 
in the disc. IBarret et alJ l|2006l ) conclude that the abrupt 
drop of the QPO coherence at about 850 Hz is due to the 
presence of the ISCO around the neutron star in this and 
other systems. 

Comparing our resu lts with Figure 8 in 
IBarret. Olive fc Miller] l|2006l ). we find that the varia- 
tion of the QPO frequency on timescales between 320 and 



2 Data in Figure [21 to [51 show no unusual change in the trend of 
the QPO frequency spread or the time derivative at around 850 
Hz, where IBarret et al.l J2006T) see the sudden drop of the coher- 
ence of the lower kHz QPO in 4U 1636—53. Similarly, there are 
no unusual features at frequencies corresponding to small inte- 
ger ratios of the kHz QPO frequencies in the spread and t'Qpo 
relations ( Abramowic z et al. l l2003l ; lBeIIoni et alj|2005l) . 
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Figure 5. Comparison between absolute values of the positive 
(red points) and negative frequency derivative (black points) for 
the lower kHz QPOs in 4U 1636-53 detected on timescales shorter 
than 64 s. 



1600 seconds contributes significantly to the QPO width in 
the frequency range between ~750 and ~870 Hz, whereas 
it does not contribute to the increase of the QPO width 
at higher QPO frequencies. Previous results show that 
variations on shorter timescales can have a large impact 
on th e width of QPO. E.g., lYu. van der Klis fc Jonken 
|200ll) . showed that in the neutron star system Sco X-l 
the frequency of the upper kHz QPO changes by 20 Hz on 
timescales of 1/6 Hz, which they relate to flux variation on 
these timescales. 

The trend of the QPO spread in Figure [2] shows a 
marginal increase of the spread at around 9 1 5 Hz, which 
is close to the frequency at which IBarret et all <|2005ch iden- 
tify the signature of the ISCO on the lower kHz QPO in 4U 
1636-53. If this apparent rise was real, it could be related to 
the idea that the high-frequency drop of the QPO coherence 
is a signature of the disc reaching the ISCO. 

Most of the models so far proposed assume that the 
kHz QPOs are produced in an optically-thick geometrically- 
thin accretion disc, and that the frequency of either the 
lower or the upper kHz QPO represents the Keplerian fre- 
quency at a given radius in the disc (e.g.. [Miller et "all 1 19981 ; 



IStella fe Vetrilll998ll Osherovich fc Titarchuklll999ri. Fo r in- 



stance, in the sonic-point model ( Miller et alJ 1 19981 ) the 
frequency of the upper kHz QPO corresponds to the ra- 
dius at which the gas starts spiralling inward at a super- 
sonic radial velocity, and changes of the QPO frequency 
reflect changes of mass accretion rat e (ram pressure and 
radiation drag). In other mode ls (e.g., IStella fc Vetril [l998l ; 
I Osherovich fc Titarchuk1ll999l ) the frequency of the QPOs 
also depends upon the radius of the disc where the QPOs 
are produced, although the mechanism by which this radius 
changes is not specified. In those models mass accretion rate 
is also likely responsible for changes of the radius in the disc 
where the QPO is produced. 
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4.1 Comparison with standard disc theory 

We can compare our measurements of vqpo of the lower 
kHz QPO with what is expec ted in the case of a thin 
disc (|Shakura fc Sunvae vl Il973h if vqpo is the Keplerian 
frequency very close to the inner edge of the disc. We 
will further assume that the disc is truncated by an un- 
specified mechanism. Acco rding to the standard disc model 
jShakura fc Sunvaevl 119731 ). matter in the disc moves radi- 
ally with a speed given by: 

v r = 0.98 a 4/5 M 3/w M" 1/4 R~ 1/4 /" 14/5 km/s, (1) 

with / = [1 - {R*/R) 1/2 ] 1/4 , where a = 0.1a s is the 
viscosity parameter, M is the mass accretion rate in units 
of 2.6 x 10 17 g/s, AL = M/1.8M Q , R is the radius of the 
inner edge of the disc in km, and 7?* is the neutron star 
radius in units of 14 km, respective!^. This is the speed at 
which matter very close to the inner edge of the accretion 
disc will move inwards one the disc has been truncated. 
Combining this expression with the Keplerian formula, we 
find that the derivative of the Keplerian frequency can be 
expressed as: 

\u\ « 300 a 4 / 5 M 3/10 M; 2/A v££ F{u kHz ) Hz/s, (2) 

where VkHz is the Keplerian frequency in units of 1000 Hz, 
and F(vkHz) is the factor f~ 14 ^ 5 expressed as a function 
of VkHz, and normalised by the same factor calculated for 
a frequency of 1000 Hz. Assuming a 1.8 Mq neutron star, 
the function F(vkHz) increases from ~ 0.6 to ~ 0.9 as the 
frequency increases from 650 to 930 Hz. From this equation 
it is apparent that the time derivative of the Keplerian fre- 
quency very close to the inner edge of a thin disc increases 
as the Keplerian frequency increases, contrary to what we 
observe in the case of the lower kHz QPO in 4U 1636-53. 
Not only that, but for nominal values of the parameters, the 
time derivative of the Keplerian frequency in the frequency 
range between 650 and 950 Hz, is about 4 orders of magni- 
tude larger than what we observe for the lower kHz QPO 
in 4U 1636-53. As it is apparent from eq. ((2]), this differ- 
ence cannot be accounted for by changing the mass and the 
radius of the NS. 

One way to reconcile the observations with the expec- 
tations fro m a thin disc is if the viscosity parameter a is 
~ 10~ 5 (cf. lBelloni et alj|l997f ). although this would not re- 
vert the trend of isqpo with the QPO frequency. The !>qpo 
vs. vqpo trend from our analysis could be in principle re- 
produced by eq.([2J) if the product ax M depended upon the 
accretion-disc radius. If we assume for simplicity that a is 
constant with radius, the trend in eq.© would be reversed if 
the mass accreted from the secondary flowed out of the disc 
at a rate that increases as the radius decreases. In order to 
fit the observations, in 4U 1636-53 the amount of accreted 
mass that crosses through the inner edge of the disc must be 
~ 5 % of the total accreted mass from the secondary. In other 
words, about 95 % of the mass accreted from the secondary 
should not flow onto the neutron star via de accretion disc. 

3 We assumed an average luminos ity value for 4U 1636—53 of 
about 0.1 Lej^ llFord et al.1 l2000h for a distance of 5.5 Kpc 
llvan Paradiis fc White! I1995T) , where -LecM is the Eddington lu- 
minosity for a 1.8-solar masses neutron star. 



Although this scenario cannot be completely ruled out, it 
is very unlikely because such a process should leave some 
trace in the spectral properties of the system, e.g., extra 
low-energy absorption, that has not been observed so far. 

The majority of the models of the kHz QPOs pro- 
pose that it is the upper kHz QPO t hat is Keplerian (e.g., 
iMiller et al.lll998l . IStella fc Vetrlll998h . Since the upper kHz 
QPO in 4U 16 36-53 is much broad er and weaker than the 
lower one (e.g.: lBarret et al.ll2005bl ). the upper kHz QPO is 
in general less significant than the lower one when both are 
measured over the same timescale. We therefore were unable 
to recover the frequency profile of the upper kHz QPO on 
the same timescales as for the lower kHz QPO to compare 
the time derivative of the upper kHz QPO wi th that of the 
Keple rian frequency in a thin disc. However, iBelloni et al.l 
|2005! ) showed that in 4U 1636-53 the frequency of the lower 
and upper kHz QPOs, vi and u u , respectively, follow a linear 
relation v u — 0.673 x V\ +539 Hz. If this relation still holds on 
short timescales, we can combine it with our measurements 
of the lower kHz QPO and estimate the time derivative of 
the upper kHz QPO frequency. We find that, similar to the 
case of the lower kHz QPO described above, the time deriva- 
tive of the frequency of the upper kHz QPO decreases as the 
frequency of the upper kHz QPO increases, which is again 
at variance with what is expected from the thin disc model. 

From the above we conclude that changes of the QPO 
frequency, at least for the case of 4U 1636-53, are not com- 
patible with the dynamics of the inner edge of a standard 
Q-disc. 



4.2 Standard disc theory plus radiation drag 

The positive and negative derivatives of the QPO frequency 
appear to follow the same trend as a function of the QPO 
frequency (see section [3] and Figure [5]). This suggests that 
there is a coupling between the mechanisms that drive the 
inner edge of the disc inwards and outwards. In this sub- 
section we explore the case in which the accretion disc is 
truncated at the so-called sonic-point radius, the position 
in the disc where the radial velocit y of the inflowing gas 
changes from subsonic to supersonic (|Miller et al.lll99a ). In 
this scenario, mass accretion rate through the accretion disc 
pushes the inner edge of the disc inwards in a similar way 
to the one we described in the previous subsection; but in 
this case, the radiation produced by the mass that eventu- 
ally accretes onto the neutron-star surface removes angular 
momentum from the disc. The position of the inner radius of 
the accretion disc is therefore set by the interplay between 
radiation drag and viscosity. Using general relativistic cal- 
culations of the gas dynamics and radiation transpor t in the 
inner edge of the accretion disc, IMiller et al.l (|l998l ) found 
that the sonic radius can be written approximately as: 

^™ R+5 (o£ufe) (luw ( in") (iF^) - 

(3) 

where R is the NS radius in kilometres, Mi is the mass ac- 
cretion rate through the disc, Me is the Eddington mass 
accretion rate, h is the thickness of the disc at a radial dis- 
tance r in the disc, and v r is the radial velocity of the gas 
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in the accretion disc. We can therefore write: 



with 



\R> 



i 0.05 M E 



4tt 



M, 



(GM) L/ ' A R~ 5/2 \R : 



(4) 



0.01 Me) 



Mi 



R 



10km 



h/r 
0.1 



10" 5 c 
(5) 

where Mi is the time derivative of the mass accretion rate 
through the disc. We can now in principle use equations (3) 
to (5) to calculate v BI vs. v BT , given the NS mass and ra- 
dius. For that we also need to specify h/r, v r and Mi. In 
systems in which accretion is the main source of radiation, 
the bolometric luminosity is proportional to the total mass 
accretion rate. To the extent that in LMXBs the X-ray in- 
tensi ty, Jx, is a good measure of the bolometric luminosity 
fe-g- jFord et alJl200Ch . M» tx J x . Putting all this together, 
and assuming that h/r and v r are constant, we find that 
z> sr oc Us/ 3 Ix /l\ ■ Using the observed variations of Jx with 
time we find that i> sr decreases by a factor ~ 2 as u Br in- 
creases from ~ 650 to ~ 900 Hz, which is comparable to 
what we found in 4U 1636-53 (see Figure [3]). 

However, it is well known that, while on timescales of 
hours kHz QPO frequency correlates with 7x, on timescales 
of days or longer the kHz QPO frequency can be the 
same even if Jx is a factor of ~ 2 diffe rent. This is the 
so-called parallel-t rack phenomenon (Me ndez et al.l 1 19991 . 
Ivan der Klisl I2OO1T ) . If the frequency of (one of) the kHz 
QPOs reflects the Keplerian frequency at the inner edge 
of the disc, and the position of the inner edge of the disc 
depends on Mi as in eq. ©, the parallel-track phenomenon 
means that there cannot be a one-to-one relation between 
Mi and Jx, and that the procedure that we described in the 
previous paragraph is incorrect. To overcome this problem, 
we proceeded as follows: Since the fre quencies of the kHz 
QPOs appear to follow a random walk (|Belloni et al.ll2005l ; 



iBarret fe Vaughan|Effll3) . we generated values of Mi follow- 
ing a random waUQ we calculated Mu and we compute d z> sr 
vs. u ST using eq. (3) to (5). Following [Miller et al.l (|l998l ). we 
took h/r = 0.1 and v r = 10~ 5 c, we assumed a NS with a 
radius of 13 km and a mass of 1.9 Mq, and generated val- 
ues of Mi between 0.01 and 0.03 Me- We further assumed 
that the Keplerian f requency at the s onic-point radius is 
the upper kHz QPO (|Miller et al.lll998l : but see below), and 
converted the sonic-point Keplerian frequency to that of the 
lower kHz QPO us i ng th e relation u u — 0.673 x ui + 539 Hz 
from lBelloni et all |2005), where vi and u u are, respectively, 
the frequencies of the lower and the upper kHz QPO. In 
Figure [6] we plot the observed relation !>qpo vs. z^qpo ( re d 
points; these are the same data as in Figure(3]) together with 
the relation between z> sr and u SI from this simulation (grey 
stars). From this Figure it is apparent that the simulation 
reproduces the observed data well. We note that although 
the range of Mi in our simulation is a factor of ~ 3 less 
than the range o f luminosities spanned by 4U 1636-53 (e.g., 
I Ford et "ah 2000), this may be accounted for if we adjust the 
values of h/r and v r accordingly. Finally, we note that we 



4 We restricted M; to be between 1% and 3% Me generating 
uniformly random distributed steps between — 10 — 4 and 10 — 4 
M E . 
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Figure 6. Comparison between the absolute value of the fre- 
quency derivative of the lower kHz QPO frequency in 4U 1636-53 
(red points) and simulated values of the frequency derivative of 
the lo wer kHz QPO using the sonic-point model from lMiller et aP 
( 1998) (grey stars). See section ^. 2l for more details about the sim- 
ulation. 



can reproduce the data equally well if we assume that the 
frequency at the sonic-point radius is the lower kHz QPO if, 
in this case, the neutron star has a radius of 14 km and a 
mass of 1.7 Mq. 

In summary, assuming that the frequency of one of the 
kHz QPOs is equal to the Keplerian frequency at the inner 
edge of the accretion disc, the comparison between the ob- 
servations and our simulations appears to lend support to 
the idea that the inner edge of the accretion disc is set by the 
interplay between viscosity and radiation drag. Our results, 
however, do not shed light about the mechanism that pro- 
duces the kHz QPOs, or whether it is the lower or the upper 
kHz QPO the one that reflects the Keplerian frequency at 
the inner edge of the accretion disc. It remains to be seen 
whether other mechanisms (e.g., magnetic drag) can also 
reproduce the observations. 
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